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Abstract 
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I ■ Microtubules are a major component of the cytoskeleton distinguished by highly dynamic behavior both in 



vitro and in vivo referred to as dynamic instability. We propose a general mathematical model that accounts 
for the growth, catastrophe, rescue and nucleation processes in the polymerization of microtubules from tubulin 



o 
o 

• I— I ■ 

dimers. Our model is an extension of various mathematical models developed earlier formulated in order to 
' capture and unify the various aspects of tubulin polymerization. While attempting to use a minimal number 
■ of adjustable parameters, the proposed model covers a broad range of behaviors and has predictive features 

m ■ 

', discussed in the paper. We have analyzed the range of resultant dynamical behavior of the microtubules by 



I changing each of the parameter values at a time and observing the emergence of various dynamical regimes, 
that agree well with the previously reported experimental data, observing the emergence of various dynamical 
regimes, that agree well with reported experimental behavior. 
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I. INTRODUCTION 



Microtubules are protein polymers made of a//3 tubulin heterodimers that form an essential part 
of the cytoskeleton of all eukaryotic cells. Besides giving structural stability and rigidity to a cell, 
microtubules play key roles in many physiological processes such as intracellular vesicle transport 
and chromosome separation during mitosis. An individual microtubule is a hollow cylinder of 25 nm 
diameter built usually from 13 protofilaments ^. While the stable subunits of microtubules are 
actually heterodimers composed of a and f3 monomers, we will refer to them for simplicity as monomers. 
These monomers exist in two different energetic states, namely bound to a molecule of guanosine 
triphosphate (GTP) or guanosine diphosphate (GDP), respectively. Only the GTP bound monomers 
are assembly competent meaning they are able to polymerize into microtubules. After the GTP 
monomers have been added to the growing microtubule, GTP bound to jS tubulin is rapidly hydrolyzed 
(dephosphorylated) to form a bound GDP subunit. It has been hypothesized since the early 1980s that 
the so-called GTP cap on the tip of the growing microtubule gives rise to the stability of the microtubule 
{27! ]. Once the GTP cap is lost, the microtubule will switch to a "collapsing" state, referred to as a 
catastrophe which is characterized by rapid depolymerization of the microtubule into its free subunits. 
However, there is also a possibility that at some point in time a catastrophically shrinking microtubule 
acquires a new GTP cap and thereby returns to the growing populatioUja situation that is referred to 



in the literature as a rescue event. We refer the reader to the papers 
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15 



30 



34| and the 



references therein for more information about this phenomenon known in the literature as dynamical 
instability. The process of microtubule polymerization dynamics both in vitro and in vivo has been 
exhaustively reviewed by Desai and Mitchison [l^. It should be stressed that under the conditions of 
high concentration of tubulin a completely different process has been observed, namely a transition to 
a regime with damped oscillations of the tubulin mass polymerized into microtubules. This occurs for 
free tubulin concentrations exceeding a critical value, that depends on the experimental conditions and 
which corresponds to saturable polymerization kinetics Dynamic instability of microtubules has 

been observed in vivo and in vitro 27|] and highlights both the non-equilibrium nature of the physical 
problem involved and the stochastic nature of the individual microtubule growth process. Interestingly, 
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es of microtubules show collective oscillations under suitable conditions first reported by Carlier 
I25I and further studied by Mandelkow et al. 



26| who showed that above a critical value 
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25 



of the free tubulin concentration, the amount of polymerized tubulin undergoes smooth oscillations 
which are damped out as the biochemical energy source in the form of GTP molecules is gradually 



depleted. Sept et al. [3J] modeled the assembly dynamics from a chemical reaction-kinetics standpoint 
and found good agreement with the experimental data. While this model captures the main features 
of the process, it is still very empirical and incomplete. Nonetheless, to introduce the main features 
of the underlying biochemical processes we have described them in Appendix [Bl These processes 
involve polymerization/depolymerization, nucleation, and catastrophe events. When at least one auto- 



catalytic reaction is added to the system, the dynamics change significantly. Sept et al. [3J], for 
example, considered an induced catastrophe event described by and incorporated into their model to 
reproduce oscillations observed in vitro. We will show in our paper that oscillations can be reproduced 
by including growth and nucleation in the model as the only nonlinear reactions. 

Deterministic mathematical models of polymer growth largely fall into two classes, depending on 
whether the length of the polymer is discrete or continuous. The latter approximation is based on the 
assumption that the typical length of a polymer is much larger than the gain in length by adding a 
e monomer unit. Partial differential equation models of this t ype have been used for example in 
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26 



34| for the case of microtubule dynamics and in ISj to study the dynamics of prion 



proliferation. 

II. THE MATHEMATICAL MODEL 

In this paper, we propose a mathematical model for the concomitant processes of microtubule 
growth, nucleation, catastrophic shrinking and rescue. Roughly speaking, microtubules of length x 
grow at a velocity ap{t) where a is a constant and p{t) the concentration of free GTP tubulin. The 
model will therefore have the form of a nonlinear transport equation. This is essentially also the 



approach taken in some of the earlier works in this area [J, [ij, Il5|, |2^. Our intention is to develop a 
generalized model that accounts for the wealth of observed behavior including nucleation, growth to 
saturation and synchronized oscillations, in addition to catastrophes and rescues. Moreover, we have 
explicitly included the presence of a GTP cap and a measure of the microtubule's age. While striving 
for completeness in the mathematical description we have also attempted to introduce a minimal 
number of model parameters. Our model contains 7 empirical parameters, most of which can be 
determined from experimental data. 

Let Y = {{x,y) G : x > y > 0} be the state space for microtubules with a GTP cap. For 
{x,y) G Y, let u{x,y,t) denote the population density of microtubules of total length x that have a 
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GTP cap and whose GDP domain has length y. This imphes that 



is the total number of microtubules in a reference volume. We assume that the GDP domain forms a 
connected set and that the remainder of the microtubule is the GTP domain of length x — y, we will 
return to this assumption in the discussion section below. If the microtubule has no GTP cap then it 
will undergo catastrophic depolymerization. In order to keep track of this process, let v{x,t) denote 
the density of microtubules of length x without a GTP domain, again in the sense of a concentration. 
In addition, we introduce the concentrations of free GTP monomers p and free GDP monomers q. 
Catastrophic depolymerization of microtubules results in the release of GDP monomers which are 
then biochemically converted into GTP monomers in a reaction sometimes referred to as "pumping" 
since it involves a biochemical energy input from the solution. The new GTP monomers then become 
available for further microtubule growth, rescue and nucleation. Fig. [1] provides a schematic depiction 
of the tubulin cycle. The equation for u is given by 

d d d 

— y, t) + {ap{t) - (3)—u{x, y, t) + 'l^u{x, y, t) = 0. (1) 

The new monomers are added at rate ap{t) and result in an increase of the overall lengths of the 
microtubules. The constant 7 > is the progression rate of the GDP zone (i.e. the speed of hydrolysis 
within the microtubule). The rate f] can be positive if occasionally a GTP bound monomer is lost from 
the microtubule. Notice that both factors ap — (3 and 7 have the dimension LT^^. The characteristic 
curves for equation ([T]) are given by 

dx dy 

In view of the second of these equations, the variable y can also be interpreted as the "age" of 
the microtubule, since hydrolysis is assumed to start immediately upon nucleation. The boundary 
condition on Fi = {(x, y) G ]R>o : y = 0} incorporates the nucleation of microtubules without a 
GDP domain. Let il){x) be the length distribution of such freshly nucleated microtubules without a 
GDP domain. This can be a uniform distribution on some interval [a;_,x+] or a (narrow) Gaussian 
distribution centered at some point Let L* = xip{x) dx be the average length of freshly nucleated 
microtubules. The nucleation reaction is generally assumed to be a nonlinear reaction although the 
exact number of monomers n that need to come together is a matter of discussion fl,0,Q. With 
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the rate of nucleation /i > 0, the boundary condition on Fi is 

ju{x,0,t) = ^p^{t)tp{x). (2) 

Thus Fi is part of the inflow boundary of the domain Y at all times. The boundary F2 = {{x,y) G 
]R>Q : X = y} may be part of the inflow boundary respectively of the outflow boundary, depending 
on whether the growth of the entire microtubule is faster than the progression of the GDP domain. 
Precisely, let R(t) = ap(t) — /3 — 7. If R{t) > 0, then we say that the system is in a growth phase. 
This allows microtubules without a GTP cap to be rescued. The boundary condition on F2 is then 
given by 

R{t)u{x,x,t) = Xv{x,t), if R{t) > 0, (3) 

where A > is the propensity of shrinking microtubules to be rescued (it has the dimension T~^). If 
on the other hand, R{t) < 0, then we say that the system is in a state of shrinking and F2 is part of 
the outflow boundary of the domain Y. Microtubules reaching the boundary F2 are transferred to the 
population of microtubules without a GTP cap v{x,t). It will be helpful later to write equation ([T]) in 
divergence form. Let 

ap{t) — (3 

7 



m 

then equation ([T]) can be written as 



Mi + V ■ = 0. (r) 

Microtubules without a GTP cap are shrinking at a rate 5 > (which represents the loss of length 
per unit time). This population has a source or a loss term, again depending on the sign of the function 
R. The equation for v is 



-\v{x,t), ifi?(t)>0 



dt ' dx 



If R{t) < the system is in a state of shrinking, and microtubules that lose their GTP cap obviously 
enter the population of microtubules without a GTP cap. If the system is in a state of growth then 
microtubules without a GTP cap are rescued and re-enter the class u through the boundary condition 
on F2. 

GDP-bound monomers are gained by catastrophes (there are no intermediate depolymerization 
products) and are converted to GTP-bound monomers by biochemical "pumping". The gain is pro- 
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portional to the number of microtubules without a GDP domain. Hence 

d 



~ ^ J v{x,t) dx — Kq. (5) 

The constant n > denotes the rate of the (first order) pumping reaction. We assume that there is 
always enough GTP available in the solution to ensure a constant rate of pumping, but in the future 
we may also include free GTP as a variable into the model. Recall that q and p are concentrations of 
"lengths" of microtubules stored in free GDP or GTP-bound tubulin monomers, respectively. In order 
to obtain the concentrations of molecules, one can calculate q = i~^q, where i is length gained by adding 
a single monomer. Flyvbjerg et al. Q, [l^ used a simple conversion where (. = 8nm/13 = 0.6 nm, 
since 8 nm is the length of a single tubulin heterodimer and 13 is the number of protofilaments in a 
microtubule. 

The population of free GTP-bound tubulin monomers is replenished by the conversion of GDP 
monomers while losses occur due to growth and nucleation of microtubules. Therefore, 

-fP=—{ap — P) / / u{x,y,t) dydx + nq — ^p^. (6) 
dt Jq Jq 

The last term in equation ([6]) indicates that n individual monomers combine during nucleation. The 
set of initial conditions is 

u{x,y,0)=u'>ix,y), 0) = p(0) = g(0) = 

We calculate the total length ||u||a;dj/dz of GDP and GTP-bound tubulin found in microtubules with 
a GTP cap by integrating with the weight x dy dx, 



\xdydx 



oo PX 



u{x, y, t)x dy dx. 



'0 JO 

The total length of GDP-bound tubulin in collapsing microtubules is given by a similar expression 

\\v{- , t)||xdx = / v{x,t)xdx. 
Jo 

The model ([I])-© conserves the total length of bound and free tubulin | |m| |a;diydx + | ["yl |xdx + ? +P, i-e. 

^ [\Ht)\\.^dyAx+\Ht)\\^dx + q{t)+p{t)^ = 0, (7) 

see Appendix |A] for details. It is also possible to calculate the total length of GDP-bound tubulin 
found in microtubules with a GTP cap as 

POO PX 

4ound-GDp(t) = / / u{x,y,t)ydydx, (8) 
Jo JO 



and so the complementary quantity, the total length of GTP-bound tubulin, is 

4ound~GTp(^) = / / u{x,y,t){x - y)dydx. 
Jo Jo 

III. PARAMETRIZATION AND NUMERICAL RESULTS 

The parameters of our model are the following, each given with their physical dimensionality 

• a, the growth rate of microtubules such that ap{t) has dimension LT~^, 

• (3, rate of loss of a GTP-bound monomer {LT~^), 

• 7, progression rate of the GDP zone {LT~^), 

• S, depolymerization rate of microtubules without a GTP cap {LT~^), 

• K, rate of the "pumping" reaction that converts GDP-bound monomers into GTP-bound 
monomers (T"^), 

• A, rescue propensity for microtubules undergoing a catastrophe (T~^), 

• fi, rate of nucleation ([p]~*^"~^)T~^), 

• n, the order of the nucleation reaction, 

• ip, distribution of lengths of freshly nucleated microtubules, 

• L*, average length of freshly nucleated microtubules (L). 

The dimensions of a and fi require some discussion. From equation (jS]) we see that [ap] = [/3] = LT~^ 
and therefore [a] = [p]~^LT^^. Hence, we are only able to use sources that report the growth rate 
as dependent on the concentration of free GTP-bound tubulin. Equation ([6]) implies further that 
[fA = [p]~*-"~^^T~^. We collect in Tabled a set of representative numerical values for the parameters 
that have been published in the literature. 

We have implemented our model numerically using matlab, the codes will be available from the 
authors upon request. To simulate the system of equations ([I])-([n]), we discretize Y into 500 x 500 cells 
where each cell has a dimension of 200 ram x 200 nm. We use an upwind scheme with adaptive time 
step for the partial differential equations and the explicit Euler method for the ordinary differential 
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equations 29j. Since the amount of free GTP-tubulin p{t) changes over time, so does the growth 
velocity of the microtubules and hence an adaptive choice of the time step is necessary to guarantee 
the Courant-Friedrichs-Lewy condition 29| . We continuously keep track of the total amount of tubulin 
in all its forms to guarantee that our numerical solution satisfies the conservation law ([7]). Using the 
length of a single unit as stated above, we can convert a concentration of tubulin as follows 



l/iM = 3.76-10^ 



=: C. 



For consistency and comparison we choose most of the parameters from only two experimental 



sources 



20 



371]. Walker et al. 371] provide experimental estimates of the polymerization rate a, the 



loss rate of GDP monomers (3 and the depolymerization rate 5. The values (summed on both plus and 
minus ends) are 

a = 2.5 /im min /iM"^ = 1.33 ■ 10"^^ L min (= ^) , 

P = 2.4 /xm min 6 = 50 /im min 
These authors also estimate the rescue and catastrophe rates as 10 min~^ and 0.36 min~^ at 10 /iM 
tubulin concentration, respectively. The latter are used to fine tune the other parameters of our 
model, i. e. 7, k, and A. We note that the integral on the boundary r2 in equation (lAip can serve as 
a definition of the catastrophe, respectively rescue frequency, depending on the sign of R. Therefore, 
we define the time-average rescue and catastrophe rates from our model as follows 

kci,t = --^ I I u{x,x,t) dx R{t) dt, if < 0, 



r2 



k — — 

"res — rp 



(9) 



u{x, X, t) dx R(t) dt, 



if R{t) > 0, 



'0 

where T is the total simulation time. For the nucleation reaction we assume that n = 2 and n 
5.9 ■ 10" VM"^ min"^ |2o|. Table 1]. 

As one initial distribution of microtubules with GTP cap we choose 



u°{x,y) 



cexp 



2.52 



(10) 



where the constant c is chosen such that ]]«''] ]a;dydx = 5 /iM, this is half the concentration of the total 
bound tubulin. The initial concentration of free GTP-bound tubulin is p*^ = 5 fiM. The remaining 
two initial data are chosen to be 0. 

In the first modeling scenario we set the parameters a.SK = X = fi = (no conversion of GDP 
monomers, no rescue and no nucleation), see Fig. [2l We plot the time evolution of microtubules in u 
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(total amount of tubulin bound in microtubules with GTP cap, solid red curve) and v (total amount 
of tubulin bound in microtubules without GTP cap, blue curve) pools, and GTP and GDP tubulin 
dimers in p (green curve) and q (black curve) pools, respectively. After an initial period of growth, 
microtubules enter the depolymerization phase due to low concentration of the GTP tubulin (green 
curve). This is because the GDP-GTP conversion process of monomers is turned off (k = 0). As a 
result, the total length of microtubules decreases. Due to lack of nucleation and rescue (/i = = A), 
this leads to complete depolymerization of all microtubules in the system, as expected. The dashed 
red curve in Fig. [2] represents the length of GDP-bound tubulin within microtubules with a GTP cap 
(see equation [H]). 

In agreement with experimental observation 2J], a solution showing damped oscillations can be 
found by introducing recycling of GDP tubulin, and nucleation and rescue processes. Using the 
parameters k = 1 min~^, = 5.9x lO^'^/iM"^ min~^ and A = 0.136 min~^, the result is shown in Fig. [31 
From now on we refer to this parameter set along with a = 2.5 fim min ~^ fiM~^, (3 = 2.4 /im min 
and 6 = 50 /im min ~^ as to the "standard" and make all changes with respect to these values. In 
Fig. [31 kinks in the curves indicate rescue events, that occur when microtubules without a GTP cap 
acquire a GTP cap under favorable growth conditions. Setting the rescue rate A to zero results in very 
similar curves without the discontinuities (not shown). The resulting rescue and catastrophe rates are 
kres = 6.7899 min~^ and /ccat = 0.1392 min~^ which are within experimentally observed ranges (see 
Table 1). Another interesting observation in Fig. [31 is that the GDP zone (of the entire population) 
follows quite closely the total length. This means that on average, microtubules maintain their GTP 
cap. Figure [H depicts the time evolution of the population density u for the set of parameters described 
in Fig. [3l As shown, the microtubules continue to grow and shrink in time. Another interesting initial 
distribution of microtubules with GTP cap is u^{x,y) = and = 10 fiM. In this scenario the 
influence of nucleation can be studied, see Fig. [51 As shown, microtubules are nucleating, polymerizing 
and then growing while the GDP-bound portion is also progressing. 

We have tested numerically the influence of the individual parameters on the behavior of our model. 
Using the initial datum from equation (fTOl) . we have varied one parameter at a time. In order to unify 
the picture somewhat, we arrange several growth curves in two-dimensional slices of the parameter 
space. In the following, we state variation of the relevant parameter as a multiple of the standard 
value, which is denoted by an asterisk. A decreased rate of monomer addition a < 0.3a* results in no 
growth of microtubules during 10 min (Fig. [HI left column). An increase in the loss rate of GTP-bound 
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monomers to 6/3* results in the formation of much shorter microtubules with a very short GTP cap 
(Fig. [6]). A reduced rate of GTP hydrolysis allows the formation of large GTP domains and delays 
the onset of the oscillation phase (Fig. [HI bottom row). On the other hand, large values of 7 are able 
to suppress microtubule formation completely (Fig. [HI top row). It is worth noting that in Fig. Owe 
observe a crossover from a regime dominated by low-level microtubule assembly to saturation kinetics 
followed by a polymerization "overshoot" with damped oscillations. The latter is a characteristic 



feature found experimentally by Marx et al. 25|] almost 20 years ago. 

Changes of the depolymerization rate of microtubules without a GTP cap 5 lead to interesting 
behavior. Strong oscillations and a prolonged existence of microtubules without GTP cap are seen 
(Fig. [7]) when 6 is small, and rescue events are more pronounced due to a longer survival of microtubules 
without a GTP cap. A very important role is played by the rate /t at which GDP monomers are recycled 
to GTP monomers. Even a very small value such as 0.1 suffices to induce oscillations, although their 
period is much longer (Fig. [HI bottom row). On the other hand, a large value such as 10k* results in a 
quick dampening of the oscillations (Fig. [HI top row). Finally, the parameter yU, influences the timescale 
for the nucleation process, but does not change the general behavior drastically (not shown). 

In addition to these variations of parameters we also varied one initial datum, namely the con- 
centration of free GTP tubulin . Now using the standard parameter set again, we see that a high 
concentration of = 10 fiM and vP{x^y) as in equation (ITOl) results in a longer persistence of os- 
cillations (Fig. [TOl) . Our model provides results regarding the nature of the oscillatory kinetics for 
microtubule assembly which is consistent with the original experimental observations reported by Car- 
lier et al. [7|[7], specifically the data shown in their Fig. 1 where the amplitude of oscillations and 
their persistence diminsh as the concentration of free tubulin decreases from 150 yuM to 50 /iM. The 
actual values of the tubulin concentration used in the experiment may not be directly comparable to 
our values due to the specific experimental conditions under which the measurements were made. 



IV. DISCUSSION 

The dynamical behavior of microtubules has attracted many investigators over^ 
to examine the microtubule behavior in regard to many biophysical aspects 



thepast tew aec; 
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'ew decades 



30 



34| . In this paper, we have proposed a new mathematical model that includes all processes taking 



place during microtubule polymerization/depolymerization, namely: growth, nucleation, catastrophic 
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shrinkage and rescue events. Our model contains the amounts of free tubuhn in both its energetic forms 
(GTP and GDP-bound) as dependent variables. This results in a nonlinear transport equation whose 
mathematical analysis will pose a serious challenge. Nevertheless, we think this will give valuable 
insight into the role of the GTP cap in maintaining microtubule stability. The earlier paper by 
Houchmandzadeh and Vallade 1^ contained expressions for free GTP and GDP tubulin, it did not 
contain a growth velocity for microtubules that was dependent on the amount of free GTP tubulin. 
GTP monomers continue to be lost, even beyond the point when they have become completely depleted 



19|, equation (9)]. This results in possibly negative concentrations of GTP monomers, avoided in our 



treatment. 

A crucial assumption of our model is that the microtubules consist of two separated and connected 
domains, a GTP domain where the growth occurs and a "trailing" GDP domain. Only very recently. 



a pape. by DMt.ov et ai U P-^ded experimental evidence ior the presence of a GTP cap in 

microtubules in vivo. The authors of [13| suggest further, that remnants of GTP tubulin left in the 
GDP domain play a role during rescue events and that growth resumes after such remnants have 
been exposed during a catastrophic depolymerization. This is clearly an exciting new development. 
Nevertheless, the remnants are likely to be very short (say one layer of GTP tubulin). It should be 
clarified that our model describes a large population of growing and shrinking microtubules where 
stochastic events in an individual microtubule have been averaged out. 

We should also point out that our model does not contain a diffusion term of the type div(D'Vu) on 
the right hand side of equation ([T]). Such a term has been a prominent feature of earlier models Q, 13], 
although it is not contained in other models 221 . The role of "length diffusion" in a mathematical 

nn 

model for linear polymer accretion has been investigated from a mathematical point of view 



Collet et al. 



Ill ] and Laurengot and Mischler 23| discussed convergence of the solutions of the discrete 



Becker-Doring system to solutions of the continuous Lifshitz-Slyozov equation under certain scaling 
assumptions. The "standard" version of the Lifshitz-Slyozov equation in 
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does not contain a 



diffusion term, implying that a population highly concentrated at a certain length initially does not 
disperse later on. 

We have performed numerical simulations using mainly parameters from only two experimental 



sources, namely 



20 



371 ]. Already with a few choices and variations of parameters, we are able to 



reproduce commonly seen dynamical behaviors, such as complete depolymerization in case of lacking 
recycling of GDP-monomers (see figure [2]) and damped oscillations in a growing population (figure [3]). 
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The parameters that have not yet been determined experimentally are - to the best of our knowledge - 
the rescue rate A and the pumping rate k, although the order of magnitude of n has been estimated by 



theoretical arguments [5|, l3J| . By varying these parameters in simulations we can predict their influence 
on the growth behavior and suggest experimental scenarios to look for. While outside of the scope 
of the present paper, we want to point out that microtubule polymerization and depolymerization is 
the target of many cancer chemotherapy drugs. The precise mechanism, by which some drugs (such 
as vinblastine and taxol) suppress dynamic instability is a topic for future modeling and experimental 
research. 



Some important points that we plan to address in the future are [21 1 



The GTP zone is generally believed to be short, a few helical rings or 40 nm at most [33j. This 
would imply that hydrolysis, under typical conditions, proceeds at roughly the same speed as the 
growth of the microtubules which in turn is dependent on the concentration of unpolymerized 
GTP-bound tubulin. Is hydrolysis of polymerized GTP-bound tubulin a catalyzed reaction? 

It is an open question whether the model developed in this paper can be adopted to describe 
situations corresponding to in vivo conditions, such as the presence of microtubule associated 
proteins (MAPs) during the polymerization process and the existence of discrete microtubule 
organizing centers. Some recent papers emphasize their role {g], [s^ in realistic representations of 
cellular processes. 
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APPENDIX A: BALANCES FOR TUBULIN AND MICROTUBULES 



We integrate equation ([1]) over the domain Y and apply the divergence theorem for weighted 
integrals 

Jn Jan Jn 

where n is the outer normal vector. With w{x,y) = x and n = I I on Fi and n = I | on 
r2 this gives, interchanging the order of integration and differentiation, 

^\\u{t)\\^dydx = J 'ju{x,0,t)xdx- J {-{ap{t) - p) + 'j)u{x,x,t)xdx 



{ap{t) — /3)u{x, y, t) dy dx 



Y 



R{t) /°° u{x, X, t)x dx, if R{t) < 



A v{x,t)xdx, if R{t) > 



OO pX 



+ {ap(t) — P) / / u{x,y,t) dy dx. 
Jo Jo 

Likewise, we integrate equation (jlj) with weight x dx. Again, after integration by parts, we have 



\ R{t) J^u{x,x,t)xdx, ifi?(t)<0 
v{t)\\xdx = -\ ^ r " / '"{x,t)dx. 

A /q v{x,t)xdx, if > j Jo 



dt 



Adding these two results and equations ([5]) and ([6]) yields equation (I7l). 

If the weight x is removed then we obtain the total number of microtubules with or without a GTP 
cap. For microtubules with GTP cap we obtain 

^\\u{t)\\Aydx = I 'yu{x,0,t)dx + {ap{t) - /3 - -f) [ u{x,x,t)dx 
= fip'^{t) / ip{x) dx + R{t) / u{x,x,t)dx, 



that is, such microtubules are gained through nucleation and lost or gained (depending on the sign of 
R{t)) through exchange with the population without GTP cap. For the latter we have 

d f 

-j:\\v{t)\\dx = -Sv{0,t) - R{t) / u{x,x,t)dx, 

i.e. microtubules without GTP cap are lost by complete depolymerization and gained or lost through 
exchange with the population with GTP cap. Taken together 

:^(||«WI|dydx + llt^(t)lldx) =/ip"W / ^{x)dx-6v{0,t), 



dt 



ri 



13 



that is, only nucleation and complete depolymerization change the total number of microtubules. 



APPENDIX B: BIOCHEMISTRY OF TUBULIN POLYMERIZATION 



The principal elements in the model of Sept et al. (3J] can be summarized by the equations that 
follow. For simplicity, the microtubule is considered as a linear polymer rather than an object of 13 
protofilaments. We shall denote a microtubule of n subunits by MT^. Note that in solution, the 
free tubulin subunits may be bound to either GTP or GDP at their exchangeable nucleotide site and 
we denote them as Tgtp and Tgdpi respectively. Note that only tubulin bound to GTP is able to 
polymerize. 

The reaction set involved in the process consists of addition, nucleation and catastrophic collapse 
is given by 

MTn + TcTP^MTn+U 

hTgtp MTn, and 

k 

MTn riTcDP- 

For simplicity, it is assumed here that all collapses are complete and the number of dimers, n' required 
for nucleation is an adjustable parameter not available directly from experiment. These equations can 
also be supplemented by the reactivation of tubulin, to make it assembly competent, which will occur 
when the concentration of GTP is high 

Tgdp + GTP ^ Tgtp + GDP. 

The free energy change associated with this reaction is less than the free energy change of GTP 
hydrolysis in solution and the difference is attributed to a structural change in the tubulin dimer. It 
is this conformational change which presumably makes assembly possible. 
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Tables and Figures 
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value 


reference 
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0.5 — 11.5 /xni niin ^ /.iM ^ 


[in 17 98 'i^ "^71 
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i.b — oo /im mm 


To y1 in AT OQ QO QC^ 0'7l 

^1 lUi XLi ^Ai 'J' 1 


7 


0.25 /xm min ^ 


ri ^1 





44 — 50 ^m min ^ 




Tl 


1-12 


f5 16 20 341 


fJ- 


5.9 X lO^M-^min"^ 


[20] 


K 


3 - 120 mmr^ 


[5,3^ 


k 

"TOS 


2-10 min"^ 


[31, 32 


^cat 


0.1 — 1 min~^ 


[31, 32 



TABLE 1: Experimental and/or computational estimates for parameters published in the literature, some of 
which are used in the model, /cres and /ccat are rescue and catastrophe rates that can be used to fine tune A, 
K and 7. 
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R(t)>0 R(t)<o 




FIG. 1: Schematic representation of the tubuhn cycle (top) and our model presented in this paper (bottom). 
In the case R{t) = ap{t) — /? — 7 > 0, the system is in a phase of growth (dashed arrows, left panel). If 
R{t) < 0, the system is in a phase of shrinking (dashed arrows, right panel). The recycling of free monomers 
is identical in both cases and only depicted in the left panel. 
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FIG. 2: Time evolution of microtubules for the parameter set a = 2.5 ^mmin" ^uM^ , (3 = 2.4/xmmin~ , 
7 = 5.4 /xmmin"^, 6 = 50 ^ummin"^ and /x = k = A = 0. Shown are the total amounts of tubulin in each of 
its forms, namely bound in microtubules with a GTP cap ||ii||xdydx (solid red curve), bound in microtubules 
without a GTP cap ||^^||a;dx (blue curve), and tubulin monomers bound to GTP p (green curve) and bound 
to GDP q (black curve). All quantities are in units of fimL^^ of tubulin. The total length of microtubules is 
decreasing at later times due to complete depolymerization of microtubules without a GTP cap and lack of 
nucleation. The red dashed curve represents the amount of bound GDP tubulin in microtubules with GTP 
cap. The inset plot shows the conservation of the total amount of tubulin, confirming equation ([7]). 
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FIG. 3: Same as Fig. [2] but for the parameter set /x = 5.9 x 10~^/iM^^ min^^, k = 1 min"^ and A = 
0.136 min~^. The arrows indicates a rescue events, the conservation of the total amount of tubuUn is shown 
in the inset plot. The resulting catastrophe and rescue rates are kca,t = 0.1394 min~^ and kj-cs = 6.7899 min~^ 
(see equation [9]) . The parameters a, (3, 7 and 6 used in figure [2] and /u, k and A used here are referred to as 
to the "standard" parameter set from now on. 
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Time: 3.65 min Time: 4.83 min 



FIG. 4: Time evolution of the population density u for the set of parameters used in Fig. [3j Plotted are 
equidensity contours beginning from a Gaussian profile. Note the change in direction of the center as time 
progresses due to the activation of the rescue process. 
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FIG. 5: Shown are again the total amounts of tubulin in each, of its forms (||'ti||xdydx 

solid red curve, ||f lUdx 

blue curve, p green curve and q black curve). Parameters are the same as in Fig. [3l but with initial conditions 
||i*°||xdj/dx = and p° = 10/iM. 




t (min) 

FIG. 6: An increase in the loss rate of GTP monomers (3 to 6/?* does not prevent formation of microtubules 
entirely, but reduces drastically their average length. Notice also the extremely short GTP cap. All other 
parameters are the same as in Fig. [3l 
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t (min) 

FIG. 7: A decrease of the depolymerization rate of collapsing microtubules to 0.25* results in irregular 
oscillations and longer presence of microtubules without a GTP cap. All other parameters are the same as in 
Fig. El 
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FIG. 8: The behavior at several points in parameter space when a and k are varied simultaneously. Each 
subfigure shows only the amount of tubulin in microtubules with GTP cap (solid red curves) and the length 
of the trailing GDP subdomain (dashed red curves). Parameters are varied with respect to the values a* = 
2.5 /xmmin^^ /xM~^ and k* = 1 min~^. The time axis covers 10 min in each case. 
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FIG. 9: As Fig. [8l but now a and 7 are varied simultaneously. 
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gives oscillations that persist for longer times. 
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